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ABSTRACT 

One of the candidates for Type la supernova progenitors, the recurrent nova RS Ophiuchi underwent the sixth 
recorded outburst in February 2006, and for the first time a complete light curve of supersoft X-ray has been 
obtained. It shows the much earlier emergence and longer duration of a supersoft X-ray phase than expected 
before. These characteristics can be naturally understood when a significant amount of helium layer piles up 
beneath the hydrogen burning zone during the outburst, suggesting that the white dwarf (WD) is effectively 
growing up in mass. We have estimated the WD mass in RS Oph to be 1.35 ±0.01 M Q and the growth rate 
of the WD mass to be at an average rate of about 1 x 10~ 7 M Q yr" 1 . The white dwarf will probably reach the 
critical mass for Type la explosion if the present accretion continues further for a few to several times 10 s years. 

Subject headings: binaries: close — binaries: symbiotic — novae, cataclysmic variables — stars: individual 
(RS Ophiuchi) — supernovae: general — white dwarfs 



1. INTRODUCTION 

Recurrent novae are binary star systems in which mass 
is transferred onto a white dwarf (WD) primary from a 
main-sequence or a red giant secondary. The eruption is 
well-modeled as a thermonuclear runaway (hydrogen shell 
flash), which occurs when a certain amount of mass (AM; g ) 
is accumulated on the surface of the white dwarf (e.g., 
iPrialnik & Kovetz"lll995l) . From their short recurrence peri- 
ods (from a ten to several tens of years) and very rapid optical 
declines, it is believed that their white dwarfs are very massive 
and c lose to the Chandrasekhar mass (e.g., lHachisu & Kato I 
2001). If the WD mass increases after every outburst, it will 
soon explode as a Ty p e la supernova (e.g., iNomotol 119821: 
lHachisu et all 1999alb1: lHachisu & Katoll2001l) . It is, there- 
fore, crucially important to know how close the WD mass is 
to the Chandrasekhar mass and how much mass is left on the 
white dwarf after one cycle of nova outburst. 

The WD envelope rapidly expands and blows winds (e.g., 
iKato & Hachisul fl994) as a result of hydrogen-flash, and its 
photospheric radius reaches a maximum and then gradually 
shrinks. Since the total luminosity is almost constant during 
the outburst, the photospheric temperature increases in time. 
We easily understand this from Stefan-Boltzmann's equation, 
£ph = 47r7? p hO'r p J 1 , where L p h, ^ p h, and r p h are the photospheric 
luminosity, radius, and temperature, respectively. The main 
emitting region moves from optical to ultraviolet, and then 
finally to supersoft X-ray, corresponding to from 7ph ~ 10 4 K 
through r p h ~ 10 6 K. The photosphere drastically shrinks at 
the e nd of the w ind phase and a supersoft X-ray phase starts 
(e.g.. lKatol[l999l) . We are able to constrain the WD mass and 
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its growth rate if turn-on/turnoff of a supersoft X-ray phase 
are detected, because they indicate the durations of wind mass 
loss (how much mass is ejected) and hydrogen shell burning 
without wind mass loss (how much mass is left). 

The sixth recorded outburst of RS Oph was discovered 
on 2006 February 12 by a J apanese amateur astronomer H. 
Narumi (Narumi et al. 2006) when it was shining at 4th mag- 
nitude. Figure Q] shows the optical de velopment of the out- 
burst (taken from Hachis u~et al. 1 120 06). The visual light de- 
cayed rapidly during the first week and then the decay grad- 
ually slowed down ("early decline phase"). It remained at 
about 10th magnitude from 40 to 80 days after the optical 
maximum, i.e., for about 40 days ("mid-plateau phase"). The 
final decline started about 80 days after the optical maximum 
("final decline phase"). It eventually decayed to 12th magni- 
tude ("post-outburst minimum"), which is a magnitude darker 
than that in usual quiescent phase. 

A very bright supersoft X-ray phase of RS Oph was ex - 
tensively observed by t he Swift XRT (e.g jBode et al.l l2006a). 
lOsborne et al. I(l2006al) reported the emergence of highly vari- 
able soft X-ray flux between 30 and 40 days after the optical 
maximum. Then the supersoft X-ray flux was stabilized at 
about 40 days and reached a maximum at about 50 days, fol- 
lowed by a linear decline between 60 and 80 days from ~ 20 
to ~ 100 counts s _1 for 0.2-10 keV (lOsborne et al. Il2006bl) . 
It rapidly declined at about 90 days. Thus the duration of su - 
persoft X-ray phase is about 60 days (Osbor ne et al. I l2006b). 

In this Letter, we calculate theoretical light curve models 
based on the optically thick wind theory and reproduce both 
the supersoft X-ray and optical lig ht curves. The main dif- 
ference from our previous models ( lHachisu & Katol f2006a) 
is to include the effect of heat exchange between the hydro- 
gen burning layer and the helium layer built up beneath the 
hydrogen burning zone. 
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FIG. 1 . — Observed X-ray count rates and theoretical X-ray flux for an adia- 
batic boundary model. Open circles: Supersoft X-ray count rates of the Swift 
XRT observations, which are extracted in the 0.3-0.55 keV energy range. 
Solid lines: Model X-ray fluxes forMwD = 1-35 Mq are plotted against time 
(days after the optical peak, i.e., JD 2453778.79; 2006 February 12.29 UT). 
Hydrogen content (X = 0. 1, X = 0.35, and X = 0.5) is attached to each curve. 
Model X-ray flux is calculated as a simple blackbody with the response func- 
tion of the Swift XRT between 0.3 and 0.55 keV, i.e., a band of 22.5-41.3 A . 
We shift vertically the X-ray flux curve to fit the peak of the count rate, be- 
cause we intend to approximately explain the duration of the supersoft X-ray 
peak but not the detailed behavior. So, no interstellar absorption is assumed 
in our model X-ray light curves. Filled triangl es and open squares: A visual 
light curve of the 2006 outburst (taken from Hachisu et al. 2006), i.e., V 
(open squares) and y {filled triangles) magnitudes. 

In $2] we present a complete light curve of supersoft X-ray 
for the RS Oph 2006 outburst. The numerical model of our 
theoretical light curves and their fittings with the observation 
are presented in |3] Discussion and conclusions follow in $4] 

2. X-RAY DATA 

We analyze Swift XRT observations available in the 
HEASARC 1 database following the standard procedures. Ob- 
servations were taken in PC (photon counting) and WT (win- 
dowed Timing) modes depending on the arrival count rate. 
Source photons were extracted from a circular region of 40" 
in PC mode and 40 x 15 pixels in WT mode. No background 
subtraction was performed as it represented less than 1% of 
the emission during the analyzed phases of the outburst and 
therefore its effect was negligible. In Figure Q] we summa- 
rize an average value of the supersoft X-ray count rates in 
every 2000s bin. A narrow energy band of 0.3-0.55 keV was 
adopted in order to avoid possible contamination by a tail of 
hard X-ray pho tons coming from high temperature plasma at 
th e shock (e.g.. | Bode et al.ll2006ah . 

Osborne et al. I d2006al) reported an oscillatory behavior in 
the Swift XRT count rate during 30 and 3 1 days after the opti- 
cal maximum, corresponding to the period in which the super- 
soft X-ray count rate first quickly rose (Fig. [T). The count rate 
began to rapidly decline at about 80 days, which is coincident 
wi th the final decline of the optical light curve as pointed out 
by lHachisu etaTl (120061) . 

The optical mid-plateau phase of RS Oph is first clearly 
identified by the y light curve in the 2006 outburst 
dHachisu et al. Il2006l) . Such mid-plateau phases are also ob- 
served in the other recurrent novae, U Sco and CI Aql, which 

1 http://heasarc.nasa.gov/ 
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FIG. 2. — Same as those in Fig[T]but for models in which heat exchange 
between the hydrogen burning zone and helium layer is taken into account. 

is interpreted as a disk irradiated by the central hot whit e 
dwarf (e.g., lHachisu et all |2000t lHachisu & Katol 120031) . 
The irradiated disk is bright when the central white dwarf is 
as luminous as that for hydrogen shell-burning phase, and it 
becomes dark when the shell-burning ends. Therefore the fact 
that the end of a supersoft X-ray phase is coincident with the 
end of a mid-plateau phase clearly indicates that this corre- 
sponds to the termination of hydrogen shell-burning on the 
white dwarf at about 80 days after the optical maximum. 

3. MODEL OF SUPERSOFT X-RAY PHASE 

We calcula te nova light curves bas ed on the optically thick 
wind theory dKato & Hachisu 111 9941) . Our theoretical model 
is consisting of a white dwarf, a disk around the white dwarf, 
and a red giant comp anion. Irradiation eff ects of each compo- 
nent are included (lHachisu & Kato 11200 lh. 

In the previous paper (Hachisu & Kato I f2006al) . we pre- 
dicted the duration of a supersoft X-ray phase of the RS Oph 
2006 outburst. Figure[T]shows our calculated supersoft X-ray 
fluxes by using the same method as in the previous paper, i.e., 
we assumed an adiab atic condition at the bo ttom of the hy- 
drogen burning layer (Kato & Hachisu 1994), so heat did not 
flow inward but only outward. We cannot reproduce the early 
emergence (day 30-40) and long duration (60 days) of the 
supersoft X-ray phase at the same time. 

Here we introduce a new treatment by which we include the 
effect of heat exchange between the hydrogen burning zone 
and helium layer. After the optical maximum, convection de- 
scends in time. So no more processed helium is carried up- 
ward but it accumulates underneath the burning zone. The 
accumulated helium layer keeps a large amount of thermal 
energy because the temperature of burning zone is as high 
as ~ 10 8 K. The temperature of burning zone gradually de- 
creases in the later phase of the outburst. Then heat flux from 
the hot helium layer becomes important in the luminosity. 
This heat flux from the thermal reservoir makes the lifetime 
of a supersoft X-ray phase much longer. Figure [2] represents 
such new light curves for X = 0.08, 0.20, and 0.35. It is very 
clear that both the earlier emergence and longer duration of a 
supersoft X-ray phase are realized at the same time. 

Our numerical model should reproduce not only the super- 
soft X-ray light curve but also the visual light curve. Fig- 
ure [3] compares the observational visual light curve with our 
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FIG. 3. — Same as those in Fig. [5] but together with the visual light 
curves. Hydrogen content in the envelope is X = 0.12 (dashed lines) and 
X = 0.17 (solid lines). These two models predict visual light curves as shown 
in the figure. Here we assume a 0.7 Mq red giant with a radius of 35 Rq for 
the companion, an irradiated disk with a radius of 47 Rq around the white 
dwarf, a binary orbital period of 455.72 days, and a binary inclination angle 
of i = 33° I Hachisu et al. 2006). 



numerical ones. In the early decline phase, free-free emis- 
sion from the optically thin ejecta dominates t he visual light 
dHachisu et al. 1 120061: lHachisu & Kato Il2006bl) . The irradia- 
tion effects of the disk and the companion red giant domi- 
nate the optical light in the mid-plateau phase, whereas the 
WD photosphere does not contribute because it shrinks to 
0.005-0.01 Rq, much smaller than the disk. The final de- 
cline started at about 80 days, corresponding to the decay of 
hydrogen-burning. 

In the model with M WD = 1 .35 M & and X = 0. 17, optically 
thick winds stopped at f w ; nd = 40 days after the optical peak. 
At this epoch, hydrogen of 0.16 x 10~ 7 M Q still exists and can 
supply more 14 days luminosity by nuclear burning. The heat 
stored in the helium layer, ~ 5 x 10 44 ergs, can also supply 
a total luminosity of 2.8 x 10 38 ergs s -1 for another 20 days 
to the supersoft X-ray phase. Therefore these energies can 
maintain the supersoft X-ray phase for a total of 34 days. As 
a result, steady hydrogen shell-burning on the white dwarf 
ends roughly at ? ste ady = 75 days. 

Another critically important value, the mass accretion rate, 
can be estimated as follows: In our model with Mwd = 
1.35 Mq and X = 0.17 (X = 0.12), the envelope mass at the 
optical peak is AMi g ~4x 1O" 6 M . This implies that the 
average mass accretion rate onto the white dwarf is M acc ~ 
2 x 10~ 7 Mq yr" 1 during the quiescent phase between 1985 
and 2006. Here we neglect a dredge-up of core materials 
mainly because carbon and oxygen was not enriched in the 
ejecta. Among the accreted matter, the wind carries away 
about 70% (50%), i.e., AM wind ~ 2.8 x 10" 6 M Q (AM wind ~ 
2 x 10~ 6 M Q ); this is much larger than the X-ray observ ational 
indication of ~ 1 x 1O" 7 M by ISokoloski et al. I (120061) . but 
roughly consistent with the infrared observational ind ication 
of - 3 x 10" 6 M Q (iDas et al. 1120061 lLane et al. 112007b . Note 
that the wind duration is f w i nd = 40 days (32 days). The resid- 
ual, ~ 1.2 x 1O" 6 M (~2x 10" 6 M Q ), accumulates on the 
white dwarf. So the white dwarf is growing at an average 
rate of M He ~ 0.6 x 10 _7 M Q yr" 1 (M He ~ 1 X 10 _7 M Q yr" 1 ) 
during 1985-2006. 
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FIG. 4. — Same as those in Fig. [5] but for different WD masses and 
chemical compositions: Mwd = 1-37 Mq, X = 0.40 (dotted line); Mwd = 
1.33 Mq, X = 0.12 (dash-dotted line); and a non-uniform chemical com- 
position model (solid line), where we increase the hydrogen content from 
Xnuc = 0.20 to X nuc = 0.40 in the nuclear burning zone only during the plateau 
phase for Mwd = 1-35 Mq. A uniform chemical composition model of 
Mwd = 1-35 Mq with X = 0.20 (dashed line) is added for comparison. 

4. DISCUSSION AND CONCLUSIONS 

In the present work, we reproduce the duration of a super- 
soft X-ray phase by introducing heat flux from the helium 
layer u nderneath the hydrogen burning zone. In our previous 
work s dKato & Hachisu Ifl994t lKato lll999t lHachisu & Katol 
1200 U l2006allbl: lHachisu et all 12006b . we had included heat 
flux from the helium layer only during the cooling phase, for 
simplicity. When the WD mass is not so massive (Mwd ;$ 
1.0 Mq) and the hydrogen content is not so small (X > 0.3), 
this effect can be neglected because only a small part of nu- 
clear energy is absorbed into the helium layer. However, we 
found that this effect cannot be neglected when the WD mass 
is so massive as the Chandrasekhar mass (~ 1 .3 Mq) and the 
hydrogen content is very small (X < 0.3), as seen in Figures 
[Jandfi] 

We have assumed that the chemical composition is uniform 
throughout the envelope for simplicity. As mentioned earlier, 
a gradient of the hydrogen content is reasonable if convec- 
tion descended in time during the rising phase of the nova 
outburst. Although the hydrogen content of RS Oph ejecta 
has not been estimated in detail, a very low hydrogen con- 
tent X ~ 0.1 was reported in the late phase of the U Sco 1979 
outburst dBarlow et alJll98lh . An interestin g trend was also 
found in the latest 1999 outburst of U Sco: lliiima I d2002b re- 
ported a rather high value of X ~ 0.6 at 16 hours aft er the op- 
tical maximum and Anu pama & Dewangan I ( 12000b obtained 
a value of X ~ 0.4 at 11-12 days after the optical maximum. 
The hydrogen content, X, may not be uniform but gradually 
increase outward. 

We assume X = 0.20 but, in the plateau phase, we increase 
the hydrogen content from X nuc = 0.20 to X nuc = 0.40 at the 
hydrogen burning zone, where X lmc is the hydrogen content at 
the nuclear burning zone. As shown in FigureH] more hydro- 
gen content extends the supersoft X-ray phase . 

In the previous work (lHachisu et al. Il2006h . we have esti- 
mated the WD mass to be 1.35 ± 0.01 Mq from the optical 
light curve fitting. Here we have estimated the WD mass and 
obtained the same results, that is, 1.35 ±0.01 Mq. A more 
massive white dwarf (1.37 Mq) can be rejected because the 
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optical light curve decays too fast in the early decline phase 
and the duration of the supersoft X-ray phase is too short (see 
Fig. @). A less massive one (1.33 M ) can also be excluded 
because we must assume a very low hydrogen content to re- 
produce the early emergence of supersoft X-ray but its super- 
soft X-ray phase lasts too long and the visual decline is too 
slow as can be easily seen from Figure [4] 

The supersoft X-ray flux from our model has been calcu- 
lated from a blackbody photosphere. Although the supersoft 
X-ray flux may not be a simple blackbody, our flux is a rea- 
sonable indication of the supersoft X-ray phase, because we 
intend to approximately explain the duration of the supersoft 
X-ray peak (not the detailed behavior). 

A few groups resolved a size of near infrared emission for 
the RS Oph 2006 outburst (IMonnier et al. Il2006h lLane et al. I 
l2007tlChesneau et al. 1120071) . iMonnier et al. I (120061) reported 
that a size of ~ 3 mas is consistent with the binary size if 
the distance is as short as 0.6 kpc dHachisu & Katoll2001l) 
but it is much larger than t he binary size when th e dis- 
tance is as long as 1.6 kp c dHjel lming et al. 1986). Re- 
centlv. lHachisu et al. I d2006l) revised their value and proposed 
a range of d = 1 .3 - 1 .7 kpc. lLane et aTl d2007l) obtained that 
the angular diameter of infrared emission increased once to 
~ 4 mas at about day 20 and then decreased to ~ 2 mas at 
about day 100. Our optically thick wind model suggests that 
the radius of near infrared (2.2/j,) photosphere is abou t 1 AU , 
estimated from equation (11) of IWright & Barlowl d 19751) . 
near the optical maximum when free-free emission dominates 
infrared continuum. The radius quickly decreases because the 
wind mass-loss rate soon decreases and free-free emission be- 
comes optically thin. This maximum radius of infrared free- 
free photosphere is smaller than the binary orbit (a ~ 1.5 AU 
at the distance of 1.6 kpc), so that a size of ~ 2-4 mas in- 
frared emission may originate from circumbinary matter as 



discussed bv lMonnier et al. I d2006l) and lLane et~aT1 d2007l) . 

Thus the supersoft X-ray light curve has led us to know 
various physical parameters of the white dwarf. In this work, 
the optically thick wind is essentially important to understand 
the nova optical light curve and its supersoft X-ray dura- 
tion because, without winds, the no va duration is too long t o 
be compatible with the observation dKato & Hachisu Ill994f) . 
Moreover, the wind duration (together with X) determines the 
amount of ejected mass and processed helium mass eventu- 
ally left on the white dwarf. Now we can conclude that the 
WD mass, 1 .35 ± 0.01 M Q , is now increasing. In fact the WD 
mass is currently increasing atMn e ~ (0.5-1) x 1O~ 7 M yr . 

We may predict the future of the white dwarf. After it 
has come through many recurrent nova outbursts, the mass 
of the helium layer reaches a critical mass and a helium shel l 
flash occurs. Its strength is weak dKato & Hachisu 1 1 19991) . 
Therefore, only a small part of the helium layer will be blown 
off in the wind, and virtually all of the helium layer will be 
burnt into carbon-oxyge n and accumulates in the white dwarf 
dKato & Hachisu I [l999h . The WD mass can grow though 
many recurrent nova outbursts and helium shell flashes. If the 
companion will supply matter by the same rate as the present 
one for another hundred thousand years, we expect an Type la 
supernova explosion of the white dwarf. 
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